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Extrinsic multiferroic hybrid structures consisting of ferromagnetic and ferroelectric layers elasti- 
cally coupled to each other are promising due to their robust magnetoelectric effects even at room 
temperature. For a quantitative analysis of these magnetoelectric effects, a detailed knowledge of 
the piezoelectric and magnetoelastic behavior of both constituents as well as their mutual elastic 
coupling is mandatory. We here report on a theoretical and experimental study of the magnetic 
behavior of BaTiOs-based extrinsic multiferroic structures. An excellent agreement between molec- 
ular dynamics simulations and the experiments was found for FesoCoso/BaTiOs and Ni/BaTi03 
hybrid structures. This demonstrates that the magnetic behavior of extrinsic multiferroic hybrid 
structures can be determined by means of ab-initio calculations, allowing for the design of novel 
multiferroic hybrids. 

PACS numbers: 75.85.+t 85.80.Jm 75.80.+q 75.70.Cn 75.30.Gw 77.80.Dj 



I. INTRODUCTION 

Multiferroic materials, 1 which simultaneously possess 
at least two long-range ordering phenomena in the same 
phase, have attracted a lot of attention in recent years 
due to their rich physics and large variety of potential 
applications. 2,3 Moreover, the mutual coupling between 
the ferroic order parameters lays the foundations for a 
variety of novel phenomena as well as enhanced func- 
tionalities and improved properties in future engineered 
material systems. 4,5 In this field, so-called extrinsic mul- 
tiferroic hybrid structures, in which ferromagnetic and 
ferroelectric compounds are artificially assembled, are 
promising candidates. They enable large and robust 
cross-coupling effects at room temperature by exploiting 
the elastic coupling between the two constituents, leading 
to extrinsic magnetoelectric effects. 6-9 Such structures 
can be realized in the form of particular composites, 10 
laminate composites, 11 vertical nanostructures, 6 and hor- 
izontal hybrid structures. 12 Concentrating on horizon- 
tal hybrid structures consisting of ferromagnetic and 
ferroelectric layers, various materials were used to re- 
alize strain-mediated multiferroic horizontal hybrids in 
the past years. 13 Even industrially produced piezoelectric 
actuators 14-17 or multilayer capacitors 18 were exploited 
to achieve extrinsic magnetoelectric effects. 

In the context of magnetoelectric effects, the electric 
control of magnetism is of particular interest. This con- 
verse magnetoelectric coupling relies on mechanical de- 
formations of the ferroelectric layer caused by the con- 
verse piezoelectric effect, or by ferroelectric/ferroelastic 
domain reconfigurations. These elastic strain changes are 
transferred into the ferromagnetic thin film clamped onto 
the ferroelectric layer, modifying its magnetic properties 
due to magnetoelastic effects. 19 In total, converse magne- 
toelectric effects in multiferroic hybrids can be described 
phenomenologically by a product tensor property includ- 
ing piezoelectric and mangetoelastic effects. 20,21 Thus, 
a detailed understanding of these effects and the elas- 



tic coupling across the interface is mandatory to predict 
magnetoelectric effects in novel multiferroic hybrid struc- 
tures. 

As a first step, we here investigate magnetization 
changes caused by magnetoelastic effects in BaTi03 
(BTO) based extrinsic multiferroic hybrid structures us- 
ing both FesoCoso (FeCo) and Ni as ferromagnetic ma- 
terials. The experimental approach is similar to Lee et 
al.. 22 They reported large changes of the magnetization 
and the electric resistivity in Lao.7Sr .3Mn03/BaTi03 
(LSMO/BTO) extrinsic multiferroic hybrids, whenever 
BTO undergoes a structural phase transition. In re- 
cent years, further experiments using LSMO/BTO, 7,23 
Fe/BTO, 24-26 Fe 3 4 /BTO, 27 - 29 CoFe 2 4 /BTO, 30 and 
Sr 2 CrRe0 6 /BT0 31 ~ 34 hybrid structures were performed. 
In these experiments, magnetization jumps were ob- 
served caused by strain changes of the respective ferro- 
magnetic thin film, which are induced by the structural 
phase transitions of BTO. However, the experimental re- 
sults could hardly be explained in a quantitative man- 
ner, since ferroelectric BTO forms a multi-domain state 
upon crossing the phase transitions without any exter- 
nal field applied to the crystal. This nonuniform strain 
state is induced into the overlying ferromagnetic thin 
film and causes changes of its magnetic properties. In 
particular, since the volume fraction of the domains in 
each ferroelastic phase of BTO may not exhibit the same 
value upon cooling down or warming up the hybrid struc- 
ture, a different behavior of the overlying ferromagnetic 
thin film was observed during decreasing and increasing 
temperature. 24,30,31 This demonstrates, it is essential to 
control both the ferromagnetic and ferroelastic domain 
configuration in order to correctly describe the manip- 
ulation of the magnetization in BTO-based multiferroic 
hybrid structures. 

By using miscut BTO substrates, we show that the 
magnetization changes of multiferroic hybrid structures 
can indeed be predicted on the basis of first-principles 
effective Hamiltonian simulations. To calculate these 
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changes, we first determine the strain state of miscut 
BTO crystals as a function of temperature by means 
of molecular dynamics (MD) simulations in Section II. 
Since the overlying ferromagnetic thin film is clamped 
onto the BTO substrate, each change of the in-plane 
strain state of the BTO crystal modifies the strain state 
of the ferromagnetic thin film. Second, the strain state 
of the ferromagnetic thin film is determined under the 
assumption of an ideal strain transfer between the BTO 
crystal and the ferromagnetic thin film in Section III. 
Knowing the strain state of the ferromagnetic thin film, 
we then show in Section IV that the magnetization can 
be calculated using a phenomenological thermodynamic 
model. To compare these calculations to experimental 
results, we fabricated FeCo/BTO and Ni/BTO hybrid 
structures and investigated their magnetic properties by 
SQUID magnetometry. In Section V we show that there 
is excellent agreement between the experimental results 
and the theoretical simulations. This demonstrates that 
our approach allows us to predict the magnetic behavior 
of existing and novel ferromagnetic/ferroelectric hybrid 
structures based on first-principles calculations. 



II. BATIO3 AS THE FERROELECTRIC 
MATERIAL 
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The use of BTO as the ferroelectric constituent allows 
us to modify the magnetic properties of the overlying 
ferromagnetic thin film not only by exploiting piezoelec- 
tric effects at constant temperature but also by using 
strain effects due to changes of the in-plane lattice con- 
stant, which occur as a function of temperature at the 
natural phase transitions of BTO. Starting in the cubic 
paraelectric state (T > T c ~ 393 K), 35 BTO undergoes 
successive structural phase transitions to three different 
ferroelectric phases as the temperature is lowered. 36 Each 
transition is accompanied by a change of the crystal- 
lographic symmetry, which causes large changes of the 
strain state of the ferromagnetic film clamped onto the 
BTO substrate. 37 The first transition occurs at the ferro- 
electric Curie temperature Tc, at which the lattice sym- 
metry of BTO changes from cubic to tetragonal. In the 
tetragonal phase, the polar axis is aligned along one of 
the pseudo-cubic (100) pc directions. Since there are six 
equivalent (100) pc axes, six ferroelectric and three fer- 
roelastic domains are possible. These domains are de- 
noted as c- or et-domains depending on weather the po- 
larization is pointing out-of-plane along the [001] pc direc- 
tion or along the in-plane directions [100] pc and [010] pc , 
respectively. Within the ferroelectric state, two more 
phase transitions occur. At T ~ 278 K, the lattice sym- 
metry changes into an orthorhombic structure and at 
T ~ 183 K, it is further reduced to rhombohedral. At 
both phase transitions, a reorientation of the polar axis 
takes place. In the orthorhombic (rhombohedral) state, 
the polar axis is aligned along one of the pseudo-cubic 
(110) pc ((lll)p C ) directions. Since the parent, paraelec- 



FIG. 1. (color online) Average homogeneous strain compo- 
nents r^ E with k — 1, . . . , 6 as a function of temperature Tmd 
calculated by performing MD simulations [cf. Ref. 38]. The 
evolution of the absolute value of r; FE simulated during cool- 
ing from 420 K without any external field applied is marked 
by dashed lines. For comparison, open symbols represent ex- 
perimental data as a function of T exp [adapted from Ref. 37]. 
The strains are calculated relative to the lattice constant of 
a = 0.3948 nm derived by LDA calculations. The tempera- 
ture of the phase transitions are indicated by vertical dashed 
lines. To reduce the number of possible ferroelastic domains 
in each phase, MD simulations are performed with an electric 
field of 4MV/m applied along the z-direction and assuming a 
miscut BTO crystal (solid blue lines). The difference between 
the pseudo-cubic [001] pc -direction ([100] pc -direction) and the 
z-direction (x-direction) was assumed to be 1° (7°). 



trie structure of BTO is cubic, it is convenient to describe 
the lattice structure of BTO as a distorted cubic struc- 
ture in the whole temperature range. 

By employing the strain changes at the natural phase 
transitions of BTO, large modifications of the magne- 
tization are expected in multiferroic ferromagnet/BTO 
hybrid structures. 22 To quantitatively calculate these 
changes, we first determine the strain state of BTO sin- 
gle crystals by performing MD simulations using the 
FERAM code developed by Takeshi Nishimatsu. 39 This 
simulation is based on a first-principles effective Hamil- 
tonian constructed from local-density approximation 
(LDA) calculations and described in detail in Refs. 38 and 
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40. For the calculations, we used a 16 x 16 x 16 supercell 
with periodic boundary conditions. 41 Furthermore, the 
system was thermalized within 60 000 time steps, and 
the properties were averaged over 60 000 steps. Since 
first-principles density-functional-thcory (DFT) calcula- 
tions normally underestimate lattice constants, a neg- 
ative pressure of p = — 5.0GPa was used in the MD 
simulations. 38 Since the ferromagnetic film in our exper- 
iments is more than 1000 times thinner than the ferro- 
electric BTO substrate, we can safely assume that the 
strain state of the ferromagnet/BTO hybrid structure is 
only determined by the BTO component. 

In Fig. 1, the absolute values of the average homoge- 
neous strain components 77™ with k = 1,...,6 (in ma- 
trix notation: 771 = e lu i] 2 = e 22 , i]3 = £33, m = 2e 2 3, 
r)5 = 2e3i, r/g = 2ei2) of BTO are simulated as a function 
of temperature Tmd during cooling from 420 K without 
any external fields. The strains are calculated relative to 
the equilibrium cubic lattice structure with a lattice con- 
stant of a = 0.3948 nm derived by LDA calculations. 38 
To compare the simulation with experimental data, liter- 
ature values published by Shebanov et al. are included 
[cf. open symbols in Fig. 1(a)]. As obvious from Fig. 1, 
the MD calculations reveal the correct sequence of phase 
transitions of BTO. However, in spite of the negative 
pressure applied, the transition temperatures are still 
underestimated and the temperature scales of the MD 
calculations Tmd and the experiment T exp do not coin- 
cide. This is usually corrected by linearly rescaling the 
temperature axis in order to adjust the theoretical and 
experimental phase transition temperatures. 42 Moreover, 
the MD calculations reproduce the experimental behav- 
ior only for temperatures T cxp < 300 K. The discrepancy 
at T cxp > 300 K can mainly be attributed to the poor 
description of thermal expansion effects in the effective 
Hamiltonian as well as to the underestimation of the lat- 
tice constants in DFT calculations. 40 Nevertheless, upon 
rescaling Tmd, the elastic behavior of BTO can be well 
described by MD calculations for T cxp < 300 K. 

It is important to note that, without applying exter- 
nal fields, a multi-domain state is expected in BTO in all 
three ferroelectric/ferroelastic phases upon cooling, since 
the different ferroelectric polarization directions in each 
phase are energetically degenerated and have the same 
probability to appear. 43 As the formation of ferroelas- 
tic domains depends strongly on extrinsic effects, such 
as sample shape, structural defects and imperfections, 44 
the simulation of the detailed strain state in each phase 
is hardly possible. To resolve this issue and to allow 
an unambiguous comparison of experiment and simula- 
tion results, a full control of the ferroelastic domain con- 
figuration, which leads to a well defined strain state in 
the BTO crystal, is mandatory. 45 This can be achieved 
by using miscut BTO crystals. In this case, the coordi- 
nate system describing the ferroelastic domains in a cubic 
reference system by [100] pc , [010] pc , and [001] pc differs 
from the directions x, y, and z describing the surface of 
the BTO crystal [cf. inset of Fig. 1(b)]. By additionally 
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FIG. 2. (color online) Schematic representation of the ener- 
getically most favorable ferroelastic domains in the rhombo- 
hedral, orthorhombic, tetragonal and cubic phase under an 
electric field applied along the z-direction using miscut BTO 
crystals. The blue arrows indicate the orientation of the fer- 
roelectric polarization. The polarization rotation induced by 
the electric field is neglected for simplicity. 

applying an external electric field along the z-direction, 
one ferroelastic domain becomes energetically favorable 
in each ferroelastic phase. The resulting ferroelastic do- 
main evolution is schematically shown in Fig. 2, where 
any polarization rotations induced by the applied electric 
field 46 are neglected for simplicity. This situation can be 
simulated by performing MD calculations with an elec- 
tric field of E z — 4MV/m, which is theoretically needed 
to ensure a single domain state in the whole temperature 
range. The difference between the pseudo-cubic [001] pc 
([100] pc ) axis and the z-direction (x-direction) was as- 
sumed to be 1° (7°). Figure 1 reveals that the cubic to 
tetragonal phase transition becomes diffuse and the tran- 
sition temperature increases by about 40 K (cf. solid line 
in Fig. 1). A slight decrease of the transition temper- 
atures between the tetragonal and orthorhombic phase 
as well as between the orthorhombic and rhombohedral 
phase can also be observed. Furthermore, the symme- 
try of the orthorhombic and rhombohedral phases are 
reduced. 47 ' 48 In spite of these symmetry changes, both 
phases are consistently labeled by the parent orthorhom- 
bic and rhombohedral phases throughout this paper. As 
shown in Fig. 1, every strain component r;^ E is uniquely 
defined in the whole temperature range. Thus, with the 
knowledge of the deformation of miscut BTO crystal, the 
magnetization changes of a ferromagnetic thin film de- 
posited on top can now be calculated employing magne- 
toelastic theory. To this end, the elastic behavior of the 
ferromagnetic thin film is determined first. 



III. STRAIN STATE OF THE 
FERROMAGNETIC THIN FILM 

In so-called horizontal multiferroic hybrid structures, a 
ferromagnetic thin film is elastically coupled to a usually 
much thicker ferroelectric substrate. Neglecting surface 
effects and assuming a perfect elastic coupling between 
both constituents, the strain state of the ferromagnetic 
thin film can be considered as well-defined and homo- 
geneous throughout the entire volume. Since the fer- 
romagnetic thin film is clamped to the thick ferroelec- 
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trie substrate, its in-plane strain components 77™, ?7 FM , 
and t]q M are totally controlled by the ferroelectric sub- 
strate: r/ FM = 77™, 77™ = 77™, 77™ = V FE . Thus, only 
the remaining components 77™, 774 M , and 77™ have to 
be calculated to determine the total strain state of the 
ferromagnetic thin film. This situation is equal to the 
general case of pseudomorphic or coherent growth of epi- 
taxial thin films on crystalline substrates. 49 By using a 
phenomenological thermodynamic model the strain com- 
ponents can be calculated employing the fact that the 
out-of-plane stress components cr™, cr™, and cr™ are 
zero, since no forces are acting on the free surface of the 
ferromagnetic thin film. 14 In the early work of Pertsev 
et a/., 50 it was shown that due to mechanical bound- 
ary conditions the equilibrium thermodynamic state of 
a thin film clamped to a rigid substrate is described by 
the thermodynamic potential g F , which is given by the 
Legendre transformation of the well known Gibbs free 
energy density g FM by 50 



~FM = g FM 
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Thus, the strain components 77™ with k = 3, 4, 5 are 
derived from cr™ = d~g FM /rjl M = 0. To determine 3™, 
the magnetic energy density u FM of the ferromagnetic 
thin film, which depends on the magnetization Mi and 
the strain state 77™, is calculated. In the following we 
restrict our discussion to situations where the magnetiza- 
tion M = M s m is well-defined and the unit vector m as 
well as the saturation magnetization M s is homogeneous 
throughout the thin film. This is fulfilled for tempera- 
tures well below the magnetic Curie temperature in the 
absence of magnetic domains. In this Stoner-Wohlfarth 
type of approach, 51 the direction of the magnetization 
can be expressed by the components of m, i.e., by the 
directional cosines m^. To take into account magnetoe- 
lastic effects, u FM (rrii, 77™) can be expanded in powers of 
rrii and 77™ for small and homogeneous deformations. 52 
Thus, u = u FM (mi,7] FM ) — it contains the lowest- 
order terms of three contributions: 
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The first term 



,FM 



(2) 

(rrij) describes the magnetic 
anisotropy, depending only on the direction of the mag- 
netization rrii. The second term w™(t7 FM ) is a function 
of the strain components 77™ and thus describes the pure 
elastic energy density of the ferromagnetic thin film. The 
dependence of the elastic constants on the magnetization 
direction rrii known as morphic effect is neglected here. 53 
The third term u™ gcl (m,i, 77™) depends on the strain 
components 77™ and the direction of the magnetization 
rrii and therefore represents the interaction between the 
elastic and magnetic anisotropy energies, i.e., the first- 
order magnetoelastic energy density. Here, the linear 
coupling between the magnetization and the mechanical 
strain is left out, since piezomagnetic or even "pseudo" 



piezomagnetic effects are not expected in ferromagnetic 
materials exhibiting a homogeneous magnetization. 54 To 
account for the shape anisotropy in ferromagnetic thin 
films with finite dimensions, an additional contribution 
«dema g (™3) = (mo/^M 2 ™ 2 . is added to Eq. (2). 55 Thus, 
the thermodynamic equilibrium state of the ferromag- 
netic thin film deposited on top of a ferroelectric sub- 
strate is determined by g FM using the magnetic energy 
density u™ = it™ + «dcma g ( m 3)- Since g FM is indepen- 
dent of the choice of axes, each term has to satisfy the 
requirements of the crystal symmetry. 56 In the follow- 
ing, we concentrate on cubic polycrystalline materials, 
which show no net crystalline magnetic anisotropy. In 
this case, the remaining energy terms it FM (?7fc = 77™) 
and u™gci( m i' Vk — r / FM ) can be expressed as € 
{1,2,3}): 56 - 58 
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Here, B denotes the isotropic magnetoelastic coupling 
coefficient, which is related to the magnetostrictive strain 
A by B = — 3(cn — ci2)A/2. 56 Furthermore, cm are 
the components of the stiffness matrix of the ferromag- 
netic material with C44 = (en — ci2)/2. 59 In general, the 
magnetoelastic coupling coefficient in ferromagnetic thin 
films deviates from the bulk value due to surface effects 
and/or the influence of strain. 60,61 Since we here use bulk 
values for B, a proportionality factor \ is introduced to 
account for any deviation in the magnetoelastic coupling 
from bulk-like behavior. The remaining strain compo- 
nents of the ferromagnetic thin film (rjf. with k = 3, 4, 5) 
can now be calculated to: 62 
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To reduce the computational cost, the magnetoelastic 
terms are neglected in Eq. (4), which yields 
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This is only valid for in-plane strains 771 and 772 larger 
than 10 -3 using ferromagnetic 3d metals like Ni or 
FeCo. By using the elastic behavior of FeCo published 
in Refs. 63 and 64, the remaining non-zero strain compo- 
nent 773 of a ferromagnetic FeCo thin film deposited on 
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FIG. 3. (color online) Evolution of the non-zero strain com- 
ponents rji M (solid lines), r/™ (dashed lines), 773 M (dotted 
lines), and 77™ (dashed-dotted lines) of a FeCo thin film de- 
posited on a miscut BTO crystal as a function of temperature 
with an electric field of 4 MV/m applied along the z-direction. 
The calculations were performed on the basis of the MD sim- 
ulation shown in Fig. 1 assuming an ideal elastic coupling 
between the FeCo thin film and the BTO substrate. The de- 
position of the unstrained polycrystalline FeCo thin film was 
assumed to take place at Tmd = 260 K on ferroelastic (a) 
c-domains, (b) ai-domains, and (c) a2-domains. 



a miscut BTO crystal as a function of the temperature 
Tmd can be calculated on the basis of the MD simula- 
tions shown in Fig. 1 with an electric field E z = 4 MV/m 
applied along the z-direction. 

In our experiments, the ferromagnetic polycrystalline 
thin films are deposited on BTO crystals at room tem- 
perature, which corresponds to Tmd ~ 260 K [cf. Fig. 1]. 
At this temperature, the BTO crystal exhibits a certain 
ferroelastic multi-domain configuration on which the fer- 



FIG. 4. (color online) Temperature dependence of the nor- 
malized magnetization projection M/M B along an external 
magnetic field H with H || x (Af x /A/ S ) (solid lines) and H || y 
(M y /M a ) (dashed lines) of a FeCo/BTO multiferroic hybrid. 
The calculations were performed on the basis of the strain 
components shown in Fig. 3 using a magnetic field strength 
of fj, H = 100 mT and an electric field E z = 4 MV/m. Depo- 
sition of the FeCo thin film on ferroelastic (a) c-domains, (b) 
ai-domains, and (c) a2-domains are considered. 



romagnetic thin film is deposited in an unstrained state. 
The detailed domain structure, i.e., the volume fraction 
of the tetragonal c- and a-domains of the BTO substrate 
during the fabrication process is unknown, since no ex- 
ternal fields are applied during the deposition. Thus, the 
elastic deformation of a FeCo thin film upon cooling is 
shown in Fig. 3 assuming a deposition of the ferromag- 
netic thin film on ferroelastic c-domains [Fig. 3(a)] as well 
as on ai-domains [Fig. 3(b)] and ci2-domains [Fig. 3(c)]. 
As expected, large changes of the strain state of the FeCo 
thin film up to 0.5% are visible while crossing the ferroe- 
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lastic phase transitions of the BTO crystal from tetrag- 
onal to orthorhombic at Tmd ~ 243 K as well as from 
orthorhombic to rhombohedral at Tmd ~ 203 K. These 
strain changes in turn affect the magnetic properties due 
to magnetoelastic effects. 



IV. MAGNETOELASTIC EFFECTS IN 
FERROMAGNET/BATIO3 HYBRIDS 

As the magnetization aligns in such a way that g FM 
takes its minimum value, the magnetization direction of 
the ferromagnetic FeCo thin film in a FeCo/BTO mul- 
tifcrroic hybrid can be calculated as a function of the 
strain state, i.e., as a function of the temperature Tmd- 
Using the magnetoelastic behavior of FeCo published in 
Ref. 64 and the correction factor x — 1-0, i.e., assum- 
ing bulk-like magnetoelastic properties, the normalized 
projection of the magnetization M/M s along an external 
magnetic field H can be calculated. 

The normalized magnetization projection with H||x 
(M x /M s ) and H||y (M y /M s ) thus obtained is depicted in 
Fig. 4. A magnetic field strength of /.i H = 100 mT was 
applied for the calculations to ensure a magnetic single 
domain state. Again, deposition of the FeCo thin film on 
ferroelastic c-domains [Fig. 4(a)], ai-domains [Fig. 4(b)], 
and a 2 -domains [Fig. 4(c)] are considered. As obvious 
from Fig. 4, large modifications of the magnetic state 
occur at the natural phase transitions of BTO even at 
a magnetic field strength of /j,qH = 100 mT. Assum- 
ing the presence of only c-domains during the deposi- 
tion, a strong magnetic anisotropy is expected in the 
temperature range of the orthorhombic phase of BTO 
[cf. Fig. 4(a)], while a finite density of a-domains mainly 
causes a huge direction dependence of the magnetiza- 
tion in the tetragonal and rhombohedral phase of BTO 
[cf. Fig. 4(b) and (c)]. This demonstrates that a different 
behavior of the magnetization is expected, depending on 
the volume fraction of the different ferroelastic domains 
present during deposition. Thus, with the knowledge of 
the domain configuration during the deposition of the fer- 
romagnetic thin film, predictions of the magnetic behav- 
ior of ferromagnetic/BTO multiferroic hybrid structures 
based on ab-initio calculations are possible. 



V. EXPERIMENT VERSUS SIMULATION 

To demonstrate the validity of the simulations shown 
in Fig. 4, the magnetic properties of ferromagnet/BTO 
hybrid structures were investigated experimentally us- 
ing FeCo and Ni as ferromagnetic materials. These hy- 
brids were fabricated at room temperature by depositing 
polycrystalline FeCo or Ni thin films with a thickness of 
50 nm onto 0.5 mm thick (OOl)-oriented BTO crystals by 
means of electron beam evaporation at a base pressure 
of 1.0 x 10~ 8 mbar. To prevent oxidation of these thin 
films, a 10 nm thick Au film was deposited in situ on top 
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FIG. 5. (color online) SQUID magnetometry measurements 
of a 50 nm thick FeCo thin film deposited on a miscut BTO 
crystal, performed under different in-plane orientations cj> of 
the external magnetic field. The measurements were carried 
out for a magnetic field strength of 30 mT (black symbols) 
and 50 mT (red symbols) at 300 K with an electric field of 
400kV/m applied along the z-direction of the FeCo/BTO 
hybrid. After magnetizing the FeCo thin film along the x- 
direction (<f> = 0°) with a magnetic field of IT, the angu- 
lar sweeps were recorded three times in both positive (full 
symbols) and negative (open symbols) direction of rotation. 
To simulate the angular dependence of the magnetization by 
means of magnetoelastic theory (solid and dashed lines), an 
energy barrier of AE/M S — 7.5 mT was assumed to account 
for incoherent switching effects of the magnetization. Fur- 
thermore, a saturation magnetization of M s = 1320 kA/m 
was used. 



of the ferromagnetic layers. Furthermore, a Au bottom 
electrode was sputtered on the backside of the BTO crys- 
tal. This enables us to apply an electric field across the 
BTO substrate along the z-direction. After the fabrica- 
tion process, the ferromagnetic/BTO hybrid structures 
were heated to 450 K, well above the ferroelectric tran- 
sition temperature of BTO, and slowly cooled down to 
room temperature under an electric field of 400kV/m, 
which results in a preferential formation of ferroelastic 
c-domains in the tetragonal phase of BTO. In the follow- 
ing, we concentrate on two hybrids using either a FeCo 
or a Ni thin film as the ferromagnetic layer. 

The so far unknown parameters, i.e., the volume 
fractions of the ferroelastic domains during deposition 
(x£ ep , cCqj P , and x^ p ) and the proportionality factor 
X [cf. Eq. (3)] can be indirectly determined by mea- 
suring the magnetic behavior of the ferromagnetic thin 
film at 300 K under different in-plane orientations of 
the external magnetic field while applying an electric 
field of 400kV/m along the z-direction of the BTO sub- 
strate, i.e., at a constant strain state of the hybrid struc- 
ture. As shown in a previous work, 65 the application of 
400kV/m ensures a ferroelastic single c-domain state in 
poled BTO crystals at 300 K. Thus our MD simulations 
largely overestimate the critical electric field required to 
obtain a ferroelectric single-domain state. This can be 
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FIG. 6. (color online) [(a), (b)] Temperature dependent magnetization M(T) measured by SQUID magnetometry for different 
values of the external magnetic field applied along the (a) x- and (b) y-direction of the FeCo/BTO hybrid structure. An 
electric field of 400kV/m was applied along the z-direction during each measurement, [(c), (d)] Simulation of the M(T) 
behavior assuming the ferroelastic domain sequence as depicted in Fig. 2. The volume fraction of the ferroelastic c- and 
a-domains during the deposition as well as the proportionality factor x was se t to x^ cp — (48 ± 5)%, = (12 ± 5)%, 

x< a^ = (40 ± 2)%, and \ = 0.42 ± 0.02 as determined by the angular dependence of the magnetization shown in Fig. 5. The 
temperatures of the experiment and the simulation are denoted by T exp and Tmd , respectively. 



explained by the inhomogeneous nucleation of domains, 
which are often pinned at defects. 48 At an electric field 
strength of 400kV/m, parts of the polycrystalline fer- 
romagnetic thin film fabricated on top of ferroelastic a- 
domains are strained, since these a-domains are trans- 
formed to c-domains at E = 400kV/m. Due to magne- 
toelastic effects, these strained regions of the ferromag- 
netic thin film cause an angular dependence of the mag- 
netization, which is different for regions of the ferromag- 
netic thin film deposited on top of o, a\- and a2-domains. 
By using SQUID magnetometry, the angular dependence 
of the magnetization projection M(cf>) along an external 
magnetic field H of a ferromagnetic/BTO hybrid struc- 
ture can be measured. 

As an example, Fig. 5 shows the angular dependence 
of M of a FeCo/BTO hybrid structure recorded at exter- 
nal magnetic field strengths of 30 mT and 50 mT. Before 
each measurement, the FeCo thin film was magnetized 
along the x-direction ((f> = 0°) using an external mag- 
netic field of 1 T. The angular sweeps were then recorded 
three times in both positive and negative directions of 
rotation. The third back (open symbols) and forth scans 
(full symbols) are displayed in Fig. 5. The (^-angles at 
which M exhibits a minimum value mark the direction 



of a magnetically hard direction. At these positions the 
measurements carried out on rotating the sample in pos- 
itive and negative direction do not coincide. The differ- 
ence as a function of the rotation direction can be at- 
tributed to an energy barrier AE describing the energy 
required to nucleate and unpin domains. 66 By employing 
the results of our MD simulations at Tmd = 260 K, which 
corresponds to T cxp ~ 300 K, performed under an elec- 
tric field of 4MV/m and OMV/m [cf. Fig. 1], M{(j)) can 
be calculated using the thermodynamic potential <7 FM . 
The energy barrier of AE/M S = 7.5 mT is chosen such 
that the angles at which the incoherent switching takes 
place are reproduced by the simulation. Furthermore, 
the experimentally determined saturation magnetization 
M s = 1320 kA/m was used. The best fit between exper- 
iment and simulation was obtained assuming a concen- 
tration of a-domains of x^ p = (12 ± 5)% and = 
(40 ±2)% during the deposition. We note that the dis- 
continuity observed in the range 130° < 4> < 210° gives 
evidence for the presence of both types of a-domains. 
Furthermore, the proportionality factor \i was found to 
be x = 0.42 ± 0.02. Using these values, the FeCo/BTO 
hybrid structure is entirely described in terms of magne- 
toelastic effects and the temperature dependence of the 
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FIG. 7. (color online) Magnetization versus temperature for a Ni/BTO multiferroic hybrid structure, (a) The measurements 
were carried out with an electric field of E z — 400kV/m and an external magnetic field H || y with different magnetic field 
strengths, (b) The simulation is carried out using a saturation magnetization of 396 kA/m, which is determined experimentally 
by measuring magnetic hysteresis loops. The volume fraction of ferroelastic c- and a-domains during deposition was set to 
xfP = 13% and = 24%. A X value of 0.80 was used. 



magnetization can be calculated on the basis of the sim- 
ulations shown in Fig. 4. 

The simulation and the experimental results of SQUID 
magnetometry measurements of the temperature depen- 
dent magnetization for different values of the external 
magnetic field applied along the x- and y-direction of 
the FeCo/BTO hybrid structure are depicted in Fig. 6. 
The data were recorded while decreasing the tempera- 
ture from 300 K to 180 K after aligning the magnetization 
into a well-defined state using an external magnetic field 
of IT. Again, an electric field of 400kV/m was applied 
along the z-direction to ensure the ferroelastic domain 
sequence in the miscut BTO crystal as shown in Fig. 2. 
In contrast to Sahoo and coworkers, 24 we do not observe 
any difference in the magnetic behavior upon cooling or 
heating the hybrid sample. This confirms that the fer- 
roelastic domain control of BTO induces a homogeneous 
strain state in the ferromagnetic FeCo thin film. Figure 6 
reveals that the simulations of M [cf. Fig. 6(c), (d)] re- 
produce the experimental results [cf. Fig. 6(a), (b)] fairly 
well. As discussed in the context of Fig. 4, the magnetic 
anisotropy visible throughout the investigated tempera- 
ture range can be attributed to the presence of ferroe- 
lastic c- and a-domains during the deposition of FeCo. 
While the magnetic anisotropy observable in the tetrag- 
onal phase of BTO is caused by ferroelastic a-domains, 
the difference of M x and M y in the orthorhombic phase 
can mainly be attributed to a finite volume fraction of fer- 
roelastic c-domains. Thus, by controlling the ferroelastic 
state of the BTO crystal during the fabrication process, 
the desired magnetic anisotropy can be adjusted in the 
whole temperature range. 

Despite of the large magnetoelastic coupling constants 
of FeCo, the changes of the magnetization projections 
M x and M y do not exceed 25%. This is mainly caused 
by the low \ value of 0.42 ± 0.02. However, \ values up 
to 0.80 were measured in Ni/BTO hybrid samples. As an 



example, SQUID magnetometry measurements and the 
corresponding simulations of the temperature dependent 
magnetization are shown in Fig. 7. The external mag- 
netic field was applied along the y-direction of the hybrid 
structure. Moreover, the elastic and magnetoelastic con- 
stants published in Refs. 67 and 68 were used for the sim- 
ulations. As obvious from Fig. 7, an excellent agreement 
between experiment and simulation was obtained for 
= (13±5)%, x*f = (24±5)%, andx = 0.80±0.05. 
In this sample, variations of the magnetization of up to 
70% can be observed at the orthorhombic-rhombohedral 
phase transition. Unfortunately, these large changes are 
not persistent and can only be achieved shortly after the 
deposition. One explanation is the presence of mechani- 
cal fatigue, which is expected to lead to two effects. First, 
mechanical fatigue can occur in the BTO crystal, which 
in turn reduces the converse piezoelectric strain. Sec- 
ond, the elastic coupling between the ferromagnetic thin 
film and the BTO crystal might be reduced. In contrast 
to ferroelectric fatigue effects, the detailed understand- 
ing of mechanical fatigue is still lacking. To clarify this 
issue is of particular technological importance for strain- 
mediated multiferroic hybrid structures. 8 However, the 
excellent agreement between experiment and simulation 
clearly demonstrates that the magnetoelastic behavior in 
ferromagnetic/BTO multiferroic hybrid structures is well 
understood. This paves the way for engineering novel 
multiferroic composite hybrid structures. 



VI. CONCLUSION 

We have investigated magnetoelastic effects in BTO- 
based multiferroic hybrid structures using FeCo and Ni 
as ferromagnetic thin film. As BTO exhibits different 
ferroelectric/ferroelastic phases, its strain state can be 
modified by simply controlling the temperature of the 
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BTO crystal. This allowed us to study one and the 
same ferromagnetic thin him under different elastic con- 
straints. Our detailed analysis showed that a ferroelec- 
tric/ferroelastic domain control is mandatory to ensure a 
predictable, homogeneous strain state of the BTO crys- 
tal. To this end, we used miscut BTO substrates and 
applied a sufficiently large electric held along the out-of- 
plane direction so that a single ferroelastic domain be- 
comes energetically favorable in each ferroelastic phase. 
This enabled us to calculate the strain state of the BTO 
crystal as a function of temperature by means of MD 
simulations based on a first-principles effective Hamil- 
tonian. Since the ferromagnetic thin film deposited on 
top of the BTO crystal is elastically clamped to the 
BTO substrate in our multiferroic hybrid structures, each 
modification of the in-plane strain components leads to 
changes of the strain state of the ferromagnetic thin film. 
Thus, with the knowledge of the elastic behavior of BTO, 
we could calculate the strain state of the overlying fer- 
romagnetic thin film using a phenomenological thermo- 
dynamic model. The calculation revealed, as expected, 
large strain changes around the natural phase transitions 



of BTO. Employing magnetoelastic theory, we were able 
to simulate the magnetic behavior of the ferromagnetic 
thin film on the basis of the elastic behavior derived from 
MD simulations. By comparing the results of these sim- 
ulations to experimental data, we clearly demonstrated 
that the magnetic properties of BTO-based multiferroic 
hybrid structures can be theoretically modeled on the ba- 
sis of first-principles MD simulations. This opens the way 
to design novel multiferroic composite hybrid structures. 
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